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Abstract 

We consider the effects of departures from minimal flavour violations (MFV) in the con- 
text of CMSSM-like theories. Second and third generation off-diagonal elements in the 
Yukawa, sfermion, and trilinear mass matrices are taken to be non-zero at the GUT scale. 
These are run down together with MSSM parameters to the electroweak scale. We apply 
constraints from fermion masses and CKM matrix elements to limit the range of the new 
free parameters of the model. We determine the effect of the departure from MFV on 
the branching ratio of b — > s 7. We find that only when the expansion parameter in the 
down-squark sector is relatively large there is a noticeable effect, which tends to relax the 
lower limit from b — > s 7 on the universal gaugino mass. We also find that the expansion 
parameter associated with the slepton sector needs to be smaller than the corresponding 
parameter in the down-squark sector in order to be compliant with the bound imposed 
by the branching ratio of r — > ^7. 



1 Introduction 



Flavour symmetries (FS), PQ-jS], are constructed using experimental data from fermion 
masses and mixings, with or without imposing a supersymmetric framework. In the- 
ories with supersymmetric FS, there is no general principle which forbids off-diagonal 
terms in the soft squared mass or the trilinear mass matrices at the electroweak (EW) 
scale. Indeed supersymmetric FS are constructed so as to be independent of the spe- 
cific mechanism for breaking supersymmetry. In addition, it is well known that in the 
minimal supersymmetric standard model (MSSM), dangerous supersymmetric contribu- 
tions to flavour changing processes can be avoided assuming flavour diagonal sfermion 
soft squared mass and trilinear mass matrices at the EW scale. However, there is no 
intrinsic reason why off-diagonal terms should be rotated away using the same unitary 
matrices which diagonalize the Yukawa matrices. In what follows, FS will always refer to 
supersymmetric flavour symmetries. 

A supersymmetric theory can be considered to be minimal flavour violating (MFV) if 
flavour transitions occur only in the charged-current sector and such that can be entirely 
determined by the CKM angles [I]. In this sense, it is natural to assume that FS can 
be postulated assuming these MFV conditions. If FS are embedded in a supersymmetric 
grand unified theory (GUT) [lj [2], the Yukawa matrices have a well defined structure 
which contains off-diagonal terms at Mqut and which evolve to the EW scale, reproducing 
the observed fermion masses and mixings. However, every possible assumption for the 
supersymmetric flavour structure at Mqut can not be compatible with this version of 
MFV. Clearly, the running of general sfermion mass squared and trilinear mass matrices 
generate angles that cannot be absorbed by the CKM elements. As a consequence, it is 
often assumed that one can impose vanishing off-diagonal terms at the GUT scale in the 
soft squared mass and trilinear mass matrices [5]. In addition, when the trilinear terms 
are diagonal, the off-diagonal terms in the squark mass matrices generated radiatively 
at the EW scale can be rotated to the MFV form defined in [5] with the same unitary 
matrices diagonalizing the Yukawa matrices. In FS with supergravity [21 El E] , the scalar 
fields breaking such symmetries (flavons) generate off-diagonal terms at the GUT scale in 
the supersymmetric sector that are correlated to those of the fermion sector in a model 
dependent way. The available model dependent freedom allows one to construct safe 
theories at the EW scale. 

Rather than determining the possible structure of the Yukawa matrices at the GUT 
scale (which is the basis for constructing FS) using only the fermion masses and mixings 
of the SM at the EW scale and hoping that the FCNC bounds can be satisfied, here we 
propose a bottom up approach for determining the parameters, both in the fermion sector 
and in the sfermion sector, of possible flavour symmetries taking into account (i) fermion 
masses with their corresponding supersymmetric corrections and mixings at the EW scale 
and (ii) measurements or bounds from FCNC experiments. 

We will begin by first discussing the flavour violating processes that we will use in 
our subsequent analysis in section 2. We will focus primarily on the branching ratio 
for b — > S7, given the current improvement in measuring its experimental value [9l [TO] . 
We discuss its computation in models without MFV in section 2 as well. In section 3, 
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we briefly discuss general features of the fermion mass matrices and the importance of 
considering their supersymmetric corrections in order to properly determine their profile 
at the GUT scale. We also point out that these theories naturally go beyond the MVF 
conditions at the EW scale. 

In section 4, we employ conditions based on supergravity and FS to determine the form 
of the soft squared mass and trilinear mass matrices that we consider. The approach we 
take here makes use of an underlying supergravity theory with an effective Kahler potential 
from which one can construct the potential and through minimization, calculate the soft 
squared masses and trilinear couplings. It turns out that in supersymmetric theories 
with non-Abelian FS, the Kahler potential [2] has fewer model dependent parameters 
than other sectors of the theory and hence can be effectively used to impose the size 
of the off-diagonal terms of soft squared mass and trilinear mass matrices at Mqut- 
Using the powerful tools of supergravity, one can then determine the structure in the 
supersymmetric sector at Mqut and test its consequences at the EW scale. The form of 
the Yukawa matrices and soft terms are developed in section 5. 

In Section 6 we present the results of the analysis. We first focus on explaining the 
impact of considering non-diagonal Yukawa matrices at the GUT scale on the running of 
off-diagonal soft squared masses and trilinear couplings. Even when off-diagonal terms 
are taken to be zero at the GUT scale, due to the dependence of the beta functions of 
these terms, they will not be zero at the EW scale, but we show that they are well below 
current experimental bounds. We then determine the profile of the Yukawa matrices at 
the GUT scale when considering their supersymmetric corrections at the EW scale. We 
do not include GUT-specific effects and assume that all GUT-scale physics has decoupled 
below the GUT scale. 

We then go on to determine the sensitivity of BR(b — * s 7) on the parameters used to 
define the Yukawa couplings at the GUT scale. Although FCNC bounds serve to restrict 
the parameters determining the Yukawa matrices, soft squared matrices and trilinear 
couplings we find a region where the contribution to b — > S7 is enhanced with respect 
to the constrained minimal supersymmetric standard model (CMSSM) with MFV and 
brings its value closer to the experimental value, thus relaxing the lower bound imposed 
on the universal gaugino mass at large tan (3. We show that this is possible while at the 
same time keeping flavour violating parameters in the lepton sector below the bounds 
imposed by BR(r — > yrf). Finally, we mention the sensitivity of the determination of the 
MSSM spectra on these new parameters. 



2 Supersymmetric flavour and CP violation prob- 
lems in FS 

It is well known that the supersymmetric soft Lagrangian, £^| SM 



MSSM 
soft 
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-QM 2 qQ ] - LM~p - uM?rf - dMj$ - eMf 

-M 2 Hu H* u H u - M 2 H H*H d - (BfiH u H d + c.c.) (1) 

introduces many new parameters that are not present in the SM. There are 106 mass 
terms, phases and mixings in the Lagrangian of the minimal version of the supersymmetric 
SM (MSSM), which cannot be rotated away by redefining the phases and flavour basis 
for the quark and lepton supermultiplets, and which have no counterpart in the SM. 

However, most of the parameters involve flavour mixing or CP violation of the type 
which is severely restricted by experimental data. In the slepton sector, the most severe 
constraint is found when Mf is not diagonal leading to the decay fi — » e 7 which could 
take place, for example, via a one-loop diagram involving a virtual bino and a wino. The 
current bound on BR(/i — > erf) < 10 -11 [12] constrains all of the parameters involved in 
the decay and in particular restricts specific terms in the squared masses of the slepton 
sector. The current bound on BR(r — > fij) ( < 7 x 1CT 8 ) [12] can also impose important 
constraints in the slepton sector. 

Another example comes from K «-> K mixing where the effective Hamiltonian gets an 
important contribution from M| 12 via a box diagram involving gluinos. Winos and binos 
can also mediate these kind of diagrams. If the squark and gaugino masses are order 
1 TeV or less, it is possible to put bounds on the parameters Am K and e K appearing 
in the neutral kaon system effective Hamiltonian which restricts the amount of down- 
strange squark mixing and CP violating complex phases that one can tolerate among the 
soft parameters. Constraints related to the second and third generations come from the 
D°, D and B°, B neutral meson systems, and the decay b — > S7. After the Higgs fields 
get vevs, the a u , a d and a e matrices defined in eq. ([1]) contribute to off-diagonal squark 
and slepton mass squared terms and hence their form is also constrained by FCNC limits. 
There are other significant constraints on CP-violating phases from limits on the electric 
dipole moments of the neutron and electron and Hg [13J. 

All dangerous FCNC and CP violation effects in the MSSM can be avoided if one 
assumes that in the soft MSSM Lagrangian, Eq. ([1]), all the mass matrices are diagonal. 
If the sfermion mass matrices are diagonal, 

M| = m|l, Ml = mil, Mj = m|l, Mf = m|l, Mf = m\\ (2) 

then all squark and slepton mixings vanish. Squarks and sleptons with the same elec- 
troweak quantum numbers are degenerate in mass and can be rotated into each other. 
There are also supersymmetric contributions to FCNC modulo mixing due to a u , ad and 
a e . Hence a further assumption usually implemented is that the trilinear terms are pro- 
portional to the Yukawa matrices 

&u = A-vXui O'd = -A-dYdi 0-e = -^-eYei (3) 

which ensures that only the squarks and sleptons of the third family can have large scalar 
cubic couplings. 

Family symmetries attempt to reproduce mixings and masses in the quark and lepton 
sector usually by breaking the family symmetry spontaneously with scalar fields, flavons, 
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and coupling different powers of them to quarks and leptons to reproduce the appropriate 
size of Yukawa couplings and hence explain the different hierarchies in fermion masses. 
However once these FS are implemented, if the theory is supersymmetric then not only will 
fermion mass terms be generated but sfermion mass terms as well and these will induce 
FCNC and CP violating effects. These effects then can be used to test different family 
symmetries. Unfortunately, unless the mechanism by which supersymmetry is broken is 
known, one is forced to make more assumptions, as in the case of the constrained version 
of the MSSM (CMSSM) discussed below. 

In this paper, we will focus our attention to the constraints imposed by the rare decays 
of b — > s 7 on departures from MFV as the parameters that we consider in this analysis 
are more sensitive to this process than to other observables such as BK(B S — > 
and AMb s - As a consequence, we will need to consider a generic method for computing 
BR(b — > s 7). QCD LO effects have been computed for the running of the chromo- 
magnetic Wilson coefficients of the effective Hamiltonian describing the decay BR(b — ► 
s 7), from the EW scale to the decay scale, //& [H]. These coefficients result of 

the virtual exchange of gluinos, neutralinos, charged Higgs bosons and charginos. 

The MFV QCD NLO contributions have been computed for the exchanges of charged 
Higgs boson and charginos in 03 H Mi HZj • In [E], the MFV MSSM NLO QCD contri- 
butions have been computed. These contributions complete the NLO calculation if the 
MFV condition is imposed at the electroweak scale. In addition, some beyond leading 
order effects (BLO), when the MFV is not assumed, have been calculated in [T9] . 

Thus in a general framework, the only complete calculations of BR(6 — > s 7) are 
available at LO. However, the present value of BR(b — > s 7) in the SM at NNLO is 
calculated [TTj 

BR(6 -> s 7) = (3.15 ±0.23) x 1(T 4 (4) 

and the present experimental value estimated by the Heavy Flavour Average Group 
(HFAG) is [10] 

BR(b -f s 7) = (3.55 ± 0.24+°;^) x KT 4 . (5) 

Here, we will not need to consider general departures from MFV at the EW scale. 
Instead, we are restricting the soft terms at the GUT scale. The expansion parameter 
of the u sector, e u , is taken to be much smaller than the corresponding parameter in 
the down sector, e^. In addition, we start with the Yukawa matrix, Y u , diagonal at the 
GUT scale, and the departure from MFV will be sensitive mainly to gluinos. We provide 
the procedure that we use to obtain BR(6 — > s 7) in the framework described below in 
Appendix [A] 

3 Fermion Masses and MFV conditions 

When flavour symmetries are constructed, the quark sector is usually fixed by precision fits 
of the CKM elements. The lepton sector is similarly determined by the UMNS elements 
although only two mixing angles out of the three mixing angles and three phases (one 
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CKM-like and two Majorana) are known. Thus, the determination of the lepton mixing 
matrix is not precise as in the case of the quark sector. Hence it is more common to 
construct flavour symmetries by fixing elements of Yukawa matrices in the quark sector 
and then adjusting parameters in the lepton sector by relating these to the quark sector 
and of course using the available information from leptons. Since only the left handed 
states enter into the electroweak mixings, only the Yukawa couplings along and above (or 
below, depending on the convention) the diagonal can be determined, so it is necessary to 
make an assumption of the Yukawa couplings below (or above) the diagonal in the Yukawa 
matrix. It is natural then to classify these choices as symmetric or non-symmetric Yukawa 
matrices. For the non-symmetric matrices, however there is less precision in determining 
the elements because the number of parameters is larger than in the symmetric case. The 
form of these matrices for the symmetric case is [25J: 



Y, 
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(6) 



(7) 



with ed = 0(\ c ) and e u = 0(\ 2 C ), where A c ~ 0.23 is the Wolfenstein parameter V us = \ c . 
The coefficients of Yukawa elements are taken to be 0(1). In general, we are free to 
independently choose either Y or Y sym for the up-type and down-type sectors. As the 
non-symmetric matrices have more parameters, we will restrict our attention here to the 
symmetric cases which will be the most restrictive. 

The extension of FS to the supersymmetric sector is somewhat arbitrary, indepen- 
dently of the many complications inherent in the MSSM. The problem rests in the lack 
of an underlying theory which relates the parameters of a FS to the supersymmetric pa- 
rameters. From a bottom up approach we can tackle this problem by using the structure 
of the running of the soft parameters in the MSSM. 

The global symmetry of the gauge sector of the SM is given by U(3) 5 = SU(3)q x 
SU(3)u x SU(3)d x and is broken only by the Yukawa couplings 



Y, 



3n x 3„ 



Y, 



3n X 3, 



Y P 



x 3 e . 



(8) 



The MFV hypothesis [1] basically constrains this breaking to the fermion sector and 
imposes Eq. (T5]) and Eq. ([3]) at some low scale to ensure small FCNCs. 

Given the dependence of the beta functions of the MSSM soft terms [26] on YjYf and 
Y f Yj, their running will alter the conditions Eq. (j2J) and Eq. ([3]) and hence the MFV scale 
is set at the EW scale. Thus, within the MSSM, the MFV hypothesis implies a strong 
restriction on the scale at which FCNCs are evaluated (e.g. //&): 



Q[Y, 



x n {Y u Y^) n \ [ X J + Xi Y u Y2] Q f 
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(9) 



As a consequence we obtain the same CKM factors for the soft masses as in the SM and 
only the flavour-independent magnitude of FCNC amplitudes are modified, for example: 

A(b -> Sj) OC [(V C KM)*ts(VcKu)tb] , ^M Bd a OC [(^CKM) 4f ,(^CKM)^ S ] 2 (10) 

Of course, the motivation for the MFV conditions is the strong suppression of BSM effects 
in flavour parameters (e.g. SM CKM fits). However, one must check carefully that the 
form of QM^Q^ in fact does not violate that assumption in the sense that it could give 
a large deviation from the MFV conditions. Once the running of Yukawa matrices is 
included this can never be the case. 

If we put the MFV conditions at the GUT scale for arbitrary initial conditions, this 
can result into catastrophic FCNCs at the EW scale. Since the main goal of FS is to 
restrict the form of Yukawa matrices, Yf, we then know the form of the terms 

YjY f , Y f Yj, a\a f and a f a\ (11) 

at the GUT scale. We can restrict then the initial conditions of off-diagonal terms of 
the soft parameters by the off-diagonal elements of the terms Eq. (jlll) . The specific way 
of correlating this can be obtained with a FS governed by a non-minimal supergravity 
theory. 

The MFV conditions in supersymmetry have been extended to include terms such as 
x 2 Y^Y"J +x^YdY^Y u Y^ +XiY u Y^YdY^ in Eq. ([9]) and a similar structure to other super- 
symmetric terms [6]. In this case, gluino contributions to FCNC can compete with those 
of chargino and charged Higgs contributions. We present examples below of such cases 
when v4q is diagonal at the GUT scale. 



4 Flavour symmetries with an underlying supergrav- 
ity theory 

4.1 Trilinear terms 

The expansion parameters appearing in the Yukawa couplings, and e n , in flavour sym- 
metries arise due to scalar fields, flavons, acquiring a vev, (6), and hence breaking the 
flavour symmetry at some scale M: 
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For example in SU (3) flavour symmetries, the Yukawa matrices have the form 
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Note that we will be assuming a symmetric form for the Yukawa matrices so that Cjj = Cji 
and we drop the subscript sym. As noted earlier, there are only 3 physical CP violating 
phases and therefore there are relations between the phases (p{j in Eq. (fT3"j) . This equation 
represents the most general parameterization of the Yukawa matrices for quarks and 
charged leptons with an SU(3) FS [2j[3]. In any supersymmetric FS, the flavon fields also 
couple to the sfermions through the trilinear terms 

MvHfQtf. (14) 

The generic form of the trilinear coupling matrices in models of supersymmetry breaking 
such as supergravity mediation, gauge mediation or anomaly mediation, is of the form 

Ma = (YfU4h- (15) 

In minimal supergravity A is a constant and hence the proportionally in Eq. ([3]) is 
achieved. Once a family symmetry is considered [2] , there are additional contributions to 
the trilinear couplings given by derivatives of the Yukawa couplings with respect to the 
flavon fields [21 18], and the trilinear terms can then be written as 

(a/k = (W(4k + 4): (16) 
where the k{j are the coefficients produced when taking the derivatives with respect to 
the flavon fields times the gravitino mass, k{j oc a/jm 3 / 2 . 



4.2 Soft squared masses 



The soft squared mass terms in general scenarios of soft supersymmetry breaking remain 
diagonal, although the diagonal terms can be different from each other. However in the- 
ories with underlying supergravity and family symmetries, the off-diagonal terms involve 
powers of the flavon fields. Of course the diagonal terms are free of such contributions 
since the terms in the Lagrangian involving the soft squared masses, Eq. (pQ), need to 
respect Hermit icity. 

Once the flavour symmetry is specified, the Kahler potential can be trivially written 

as 

K = X^V^ijW; ?P = ur, d R , e R , u R , Q l , L l , (17) 



K f M) 



iil>) + d{X p ,i/,)X>Xl 
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a.3 



c(6 a ,^) + d(e a ,x p ,^)xvxl\ .(18) 



M\e a ) 

Here 9^ are the flavon fields and X p are the fields breaking supersymmetry. The coefficients 
c and d are a-priori un-related coefficients that can be different for each different flavon 
and X field. In abelian flavour symmetries, they are different for each fermion species 
but in non- Abelian symmetries the coefficients of one species to another are often related. 
For example in SU(3) flavour models, [21 El [25], we have structures like 

/ n o \ 

r 2 O ' 
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(19) 
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where it is important to stress that due to the vacuum alignment of the flavon fields, 
[21 E] there are no D terms contributing to the soft-squared masses, which often present 
dangerous contributions to FCNC [27J. In what follows, we adopt the following form for 
the soft-squared masses 
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(20) 



5 Determination of the flavour structure at the GUT 
scale 

5.1 General assumptions 

Having specified the assumptions that we make based on an underlying supergravity 
structure, we next describe how these structures are realized at the electroweak scale and 
how we might implement experimental constraints on the structure at the GUT scale. In 
a complete theory of FS, the full superpotential and Kahler potential would be specified, 
including all the coefficients appearing in Eq. ffl2l . Eq. ffl6|) and Eq. ffl8l) . One can then 
impose those initial conditions at the GUT scale and run those parameters down to the 
EW scale. Lacking a complete theory we can at best begin with the general form of the 
Yukawa matrices, Eq. ([3]), and the supergravity assumptions described in the previous 
section at the GUT scale and apply the constraints from FCNC, mainly BR(b — > S7) to 
limit the possible values of the parameters of the flavour structure for Yukawa matrices 
and soft terms at the GUT scale. To this end, we implement an iterative procedure, like 
the one followed in analyses of the constrained MSSM (CMSSM) [281129]. 

In the CMSSM, gaugino masses are assumed to be unified at the GUT scale with value 
m 1( / 2 , as are the soft scalar masses with value mo and trilinear couplings with value A . 
Typically, the flavour structure for the scalar masses and trilinear terms is assumed to 
be diagonal at both the GUT and EW scales. That is, the running of the off-diagonal 
terms is also neglected. With this set of boundary conditions at the GUT scale, we can 
use the radiative electroweak symmetry breaking conditions by specifying the ratio of the 
two Higgs vacuum expectation values, tan/3, and the mass, Mz, to predict the values of 
the Higgs mixing mass parameter, /x and the bilinear coupling, B. The sign of [i remains 
free. 

The CMSSM is subject to a number of phenomenological constraints. These include 
the LEP limits on the chargino mass: m x ± > 104 GeV [30] . on the selectron mass: 
mg > 99 GeV [31] and on the Higgs mass: > 114 GeV [32]. The former two constrain 
m 1( / 2 and tuq directly via the sparticle masses, and the latter indirectly via the sensitivity 
of radiative corrections to the Higgs mass to the sparticle masses, principally m^- b . The 
Higgs limit imposes important constraints principally on m^j particularly at low tan/3. 
Another constraint is the requirement that the branching ratio for b — > s'j is consistent 
with the experimental measurements J9j[10]. These measurements agree at the l-a level 
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with the Standard Model, and therefore provide bounds on MSSM particles [17], such as 
the chargino and charged Higgs masses, in particular. Typically, the b — > s'j constraint is 
more important for fi < 0, but it is also relevant for fi > 0, particularly when tan f3 is large. 
The constraint imposed by measurements of b — > S7 also excludes small values of mm. To 
apply the constraint from b — > 37 to the supersymmetric parameter space, we convolve the 
experimental result given in Eq. [5] with the theoretical calculated result and uncertainty. 
We add linearly the uncertainty generated from the variation of the input matching scale 
to the SM error given in EqlU We also add linearly the systematic uncertainties in the 
experimental result given in Eq. [5] which could allow central values of b — > 57 significantly 
outside the range inferred from the statistical error in Eq. [51 For example, in the CMSSM, 
the 95% CL lower limit to at high tan/5 occurs when the calculated central value of 
the branching ratio of b — ► 57 is as low as 1.9 x 10~ 4 [33]. Finally, there are regions of 
the (mi/2, m ) plane that are favoured by the BNL measurement [34J of — 2 at the 2-a 
level, corresponding to a deviation from the Standard Model calculation [35] using e + e~ 
data. 

If R-parity is conserved, the lightest supersymmetric particle becomes a cold dark 
matter candidate [36] and therefore the model becomes subject to cosmological constraints 
on the relic density of dark matter (usually the lightest neutralino) [29J. In the context 
of the ACDM model, the WMAP only results indicate [37] 

VcDMh 2 = 0.1042^080 (21) 

or a 2a range of 0.0882 - 0.1204 for QcDAih 2 . This range of densities forces one into 
relatively narrow regions in a mi/2, mo plane when tan/? and A$ are fixed. 

In this paper, we will restrict our attention to models with fixed tan/3 = 50 with 
\i > 0. This is done to enhance somewhat the importance of b — > 57 constraints. We 
take m t = 171.4 GeV and m h (m b ) = 4.25 GeV. In the CMSSM with tan/5 = 50 and 
Aq = 0, the Higgs mass constraint requires m 1( / 2 > 350 GeV and is similar to the constraint 
from b — > S7 for this choice of tan/3 and Aq. At larger values of m.1/2, the value of mo 
is constrained to lie in a relatively narrow range to ensure that the co-annihilation of 
neutralinos and staus are effective in reducing the relic density to the WMAP range. 
At very large values of m^, s-channel annihilations of neutralinos through heavy Higgs 
scalars and pseudo-scalars control the relic density and larger values of mo are allowed. 
The parameter choices used in the results presented in the next section have been made 
such that the the relic density is within the WMAP range in the CMSSM. We also note 
that to explain the g^ — 2 result at the 2-a level, places an upper bound on mi/ 2 of 
approximately 900 GeV. We have not considered here parameter choices with very large 
values of m typically associated with the focus point as this part of the plane is not very 
sensitive to FC processes. 

5.2 Form of the Yukawa matrices and soft terms 

Our goal in this paper is to test the effect of including off-diagonal terms in the scalar 
and trilinear matrices. We will restrict our attention to the mixing between the second 
and third generations only (mixing with the first generation will be of higher order in the 
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expansion parameters e u or as discussed in the previous section). At the EW scale, 
we set up off-diagonal entries of using the CKM mixings. We also assume the same 
mixings for the charged lepton sector. At the GUT scale we choose the following form of 
the Yukawa matrices, which contain 10 parameters: 




}••„ = I c u 22 e 2 u o i . t;, 



/ 
4& 4A I , Y e = I cf 2 e* cf 3 e* I . (22) 










e ,2 
22 t d 


r e 2 


e ,2 
23 t d 


°33 



In the up sector, we have the masses of the top and charm quarks to determine two of 
the three parameters in Y u , thus we have just one free parameter which we choose to be 
e u . The d and e sectors are linked because we are choosing to generate the mixing in the 
lepton sector using the Vckm matrix: 

Yd = ^CKM^^CKM' Y e = K^KM^e^CKM) (23) 

where Yf are the diagonal matrices of the bare Yukawa couplings at the EW scale. In the 
e and d sectors combined, there are seven parameters. We can use the masses of fi, r, and 
the masses of the strange and bottom quarks to eliminate four parameters. An additional 
parameter is eliminated by the Vckm matrix, leaving us with two free parameters which 
we choose to be and c 22 . Thus choosing 

e u , e d , and c d 22 (24) 

as free parameters at the GUT scale allows us to determine the remaining coefficients 

c 22> c 33' C 23' C 33' C 22> c 23' allC ^ C 33' (25) 

The fits are done at the appropriate low energy scale and run up and down iteratively to 
obtain a consistent solution given our set of input parameter values. 

Although we are choosing e u and as free parameters, their range is restricted by 
requiring convergence of the running of the RGEs. Outside this range, Yukawa couplings 
tend to diverge while running down to the EW scale. This is reminiscent of the restriction 
to tan (3 in the CMSSM. The allowed range can be estimated from the running of the mass 
ratios m s /mb and m c /m t , which are directly correlated to the values and e u respectively 
by Eq. (38). These ratios increase from the EW to the GUT scale in the MSSM. We 
use an analytic form [38] for the ratios of Yukawa couplings derived from one-loop RGEs 
for quark mixings and solved in the presence of the heaviest Yukawa generation. In this 
approach the evolution parameter, x is simply related to the evolution of the mass ratios 
rris/mb and m c /m t by 

(M G y y2 t/ 1 ^ 2 _ n o (m s /m b )(M G ) (m c /m t )(M G ) _ 3 



X VmJ ^°- 8 ' (m s /m b )(M z ) X > (m c /m t )(M z ) X ' (26) 



Then at Mqut, Q ~ \jfn s / m b ~ 0.1 and e u ~ Jrajra t ~ 0.05 and hence the range for 
and e u should be restricted to be close to 0.1 and 0.05 respectively. 
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With the above parameterization of Yukawa matrices then we can parameterize the 
soft masses as follows: 



M% = 



Ml 



a„ = 



Ml = 



( 1 
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X 23 € Su 



7 u ,2 
Z 22 e Su 
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rrir 



A), 



a P . = 



7 e 2 
z 22 t Se 



V z 2 3 4e 



Z 23 € Se 
Z 33 



Ac 



(27) 



In the spirit of the CMSSM, we can simply choose all of the x{j = 1, eg/ = e/, and set 
the trilinear couplings to be aligned with the Yukawa matrices, such that z{j = c{j. We 
will make this simplification in sections 6.1 - 6.3. 

Later, in section 6.4, we will use the generality obtained when esd, tse, and es u are 
assumed to be unrelated to and e u to set constraints on these parameters from b —>■ s 7. 
We can test departures from the CMSSM by effectively taking coefficients x{j different 
from unity, but that is clearly beyond the scope of the present study. The matrices are 
then simplified to: 
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c 23 t Sd c 33 
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A , 



(28) 



Note that for this case the trilinear couplings (~ z lj^f) wn l n °t be exactly aligned with 
the Yukawa couplings but will still be proportional to them and effectively we are taking 



c{j4d/4 f o r (hj) ^ (3,3) and z J 33 



J 

c 33- 
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6 Results 



To get a feel for the range of parameter values of the coefficients c{j and e u and e d we first 
consider a specific point in the susy parameter space with mi/2 = 520 GeV, mo = 330 
GeV, Aq = , and tan/3 = 50. We also consider a case with Aq = —mu^ and we restrict 
our attention to fi > 0. In the CMSSM , this point corresponds to one lying on the WM AP 
co-annihilation strip [2H] so the that the relic density lies within the WMAP preferred 
range. The Higgs mass is suitably large ( > 116 GeV), the susy contribution to b — ► 57 is 
small, and the contribution to (g — 2) M accounts for the discrepancy between theory and 
experiment. 

The goal of the analysis is to constrain the allowed parameter space both from flavour 
changing neutral currents and from the parameters determining the appropriate values 
of fermion masses at electroweak scale. As noted earlier, our free parameters are c\ 2 and 
(e u ,€d). We have set cf 2 = 0.7. The other c{j are then fixed by the low energy fermion 
masses and mixings. Since we expect the parameters e u and e d to describe the structure 
of the Yukawa matrices, we do not want to rely on fine tuned coefficients cf -. Indeed, 
would preclude a value of as small as 0.1 if e d is also taken of order 0.1. 

6.1 Running of off-diagonal soft parameters 

Since we work in the context of the CMSSM, it is important to first establish the degree 
of running of our off-diagonal terms between the GUT and EW scales. The 1-loop RGE's 
for (Mfi)ij, (Mj )ij and a d are given b}0 [26] 



/?« = (Mi + 2M%)Y^Y U + (M~ + 2Mjj)Y}Y d + W^Y U + YjY d 



+ 2Y^M 2 u Y u + 2Y}M 2 d Y d + 2a d a d 
P% = (2Mj + 4M 2 Hd )Y d Yj + AY d MlY^ + 2Y d Y}Ml 

d 

+ 4a d 4-y^|M 3 | 2 -^ 1 2 |M 1 | 2 + ^ 1 2 5 
/?W = a d hr(3Y d Y} + Y e Y^ + 5Y}Y d Y^Y u - £ 

[ a=l,2,3 

+ Y d \ Tr{6a d Yj + 2a e Y e t ) + 4Y}a d + 2Y^a u - 2 £ u a g 2 a M a ) (29) 
[ 0=1,2,3 

where u a = 7/9,3,16/3 are numerical coefficients, g a are the gauge couplings and M a 
are the masses of the gauginos. We have not shown (though we do include them in our 
calculations) the beta functions for M|, M 2 , and Mf. /3i7 ? is similar to /3jj7 2 with the 
interchange u <-> ti and different contributions from the gauginos. The running in the u 



1 Note that we employ a sign convention for a which differs by a relative sign from that of 
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sector is in general small due to our assumption that the expansion parameter in the up 
sector is small in comparison with that of the down sector. The running of M| (Mf ) is 
similar to that of M~ (Mj) due to their common dependence on e u (e d ). Here, we will 
describe the running of (M~)jj, (Mj and a d , but the running for (M|)y, (M~ R )ij and 
a e , has a similar behaviour. 

Y d is not diagonal at the GUT scale, however, Y u is, and although at two loops the 
evolution of Y u produces off-diagonal elements we neglect them here, due to the smallness 
of e u . Therefore, the most important contributions to the beta functions of the off-diagonal 
terms of the the soft squared masses in Eq. (|29|) are found in the terms 



(MffiYfa, (YjYfM 2 ^, (YjMjYf)^, {a\a d ) %J 



(30) 



For the off-diagonal trilinear terms, important contributions arise from the terms: 

(a d Y^Y d ) tJ , {Y d Yla d ) iv {Y d Ti{a d Y ] d )) lv (Y d Tr(a e Y e % r (31) 

Of course we need only consider the case for / = d and = (2,3). Expressing Y d in 
terms of e d , the contribution to the beta functions of the soft squared masses is 



(M 2 6 YjY d ) 
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d 32 

C 22 C 23^d 22 6 d + C 23 C 33^J 33 6 d — C 23 C 33 e d m 



r d r d /If 2 - f 2 
C 23 C 33 JW Q 33 £ rf 
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r d r d 2 2 
c 23 c 33 fc d" i 



C 22 C 33^d 23 C 23 C 33^d 33 



„d „d „2„ 



7 d 7 d f 2 A 2 







(32) 



where we have kept only like terms to lowest order in e d . The expressions on the right-hand 
side of the arrows are the remaining non-zero contributions when the off-diagonal terms 
in eqs. (T5|) and (j3J) are set to zero. From Eq. ff27|) we can see that since M~ 22 >> M? 23 , 
even when the off-diagonal terms of Eq. (|30|) are set to zero, the leading terms ~ M~ ..e d , 

of Eq. fl32|) remain, except for the term a) d a d , which vanishes. For the contribution to the 
beta functions of the trilinear couplings we have 



(a d Y}Y c 



d)23 
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d d d 4 A 
x; 22 c 23 c 33 t cZ^ 1 



{YdY d a d ) 2 3 

Y d Tr(a d Yj 



Y d Tr(a e Y?) 



23 



23 



z 23 e d Jii o{Y d Y d )22 + z d 3 A Q (Y^Y d ) 2 3 



d d d 2 a 
x; 33 c 23 c 33 t d^ 1 



d d 6 d a 
L 22 c 23 t d z 22^ 1 



O n d z c 6 d a , d d 2d A 
zc 23 fc d z 23 /1 i c 23 c 33 fc d' 5 33 /1 



d d 6 d | d d 2 d 
c 22 c 23 fc d^"22 ' c 23 c 33 fc d z 33 



A 



e d c b e a , o d e b e 4 , d e 2 e 4 
L 22 c 23 t d z 22^ 1 ~r ■ it -23 c 23 t d^23^ 1 "I" t -23 c 33 t d z 33^ 1 

r e r d 7 e -I- r d r e c 2 7 e 1 /l„ 
c 22 c 23 fc d^22 ^ c 23 c 33 fc d z 33 ^0 



(33) 



In Eq. ( 1331) we see again that when the off-diagonal terms of the trilinear couplings a djl 
are set to zero the leading terms of the original expression do not vanish. 
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In Figure [JJ, we show the running of the off-diagonal soft masses. We begin by dis- 
cussing the case for A = 0. Here, we have chosen e d = 0.1 and e u = 0.05 as well as 
c| 2 = 0.7 which will be used throughout this section. The running of the soft masses in 
this case is shown by the dashed (when the off-diagonal terms of the soft squared masses 
are set to zero at the GUT scale) and solid (when they are non-zero) blue curves. In the 
first panel of Fig. [TJ we show the running of M% as a function of the renormalization 

scale Q. When the off-diagonal terms are turned off, the running starts at M~ = 
as seen by the dashed (blue) curve. The dominant contributions to the running of are 
suppressed by e 2 d and proportional to (4 3 cf 3 . As a consequence, the running of M? 23 is 
relatively significant. When the off-diagonal elements are turned on at the GUT scale, the 
initial value for M? 23 is about (16.5 GeV) 2 (as seen more clearly in panel b for M? 23 ). 
On this scale, the difference in the running is relatively minor as all of the leading contri- 
butions to these expressions are the same. As a result the evolution of these terms does 
not change drastically according to the different initial conditions at the GUT scale. 

Also shown in Fig. [TJi is the analogous running when A 7^ 0. Here we show the 
running for the specific choice A = —m\n = —520 GeV. Note that for this case, we have 
adjusted mo to 370 GeV in order to maintain the correct relic density. When A 7^ 0, 
the evolution of soft squared masses will have contributions from the terms aW/ when 
the off-diagonal terms start with a non-zero value at the GUT scale. As the dominant 
contribution to the running of Ml 23 comes from the gaugino (mostly gluino) terms in 

Eq. ( f29|) . a^cif tends to lower (slightly) the running as seen by the dash-dotted (when 
off-diagonal terms are turned off at the GUT scale) and dotted (when they are turned 
on) (red) curves. 

In contrast, in panel b) of Fig. [TJ we show the running of Mj 23 . The beta function 
for Mj 23 is determined by Y u or a u rather than Y d and as a consequence the running is 
minimal. As expected from the beta functions, the running of Mj 23 is also significant as 
seen in panel c). We also show the running of M| 23 in panel d). As noted earlier, the 
beta function for M| 23 is similar to that of M? 23 though the running is smaller due to 
the lack of gluino contributions. 

Although the running of (MD23 is the most significant of those shown in Fig. [TJ in the 
computation of BR(b — > s 7), discussed below, this effect is sub-dominant. The reason 
is that the contributions from the charged Higgs boson and the charginos is much more 
important than the contribution from gluinos. We can see from Eq. (1791) that only the 
diagonalization matrix of Mj, K u , is relevant for the computation of the magnetic Wilson 
coefficients where the charginos and charged Higgs boson are the virtual particles. The 
diagonalization matrix of Mj, K d , is relevant only for the gluino contributions. This is 
because they enter into the computation of chromo-magnetic Wilson coefficients arising 
from the virtual exchange of gluinos [T4] . 

For A = the renormalization group equations generate a contribution for a d of order 
—Yd J2 a u aM a at one loop. In addition, when A Q = 0, the running of a d is independent of 
whether off-diagonal terms of soft parameters are set to zero at the GUT scale. In Figure 
EJ we show the running of 6^23 in panel a) and a e2 3 in panel b). The solid (blue) line 
corresponds to the case where A = 0. In both cases the running is very slight. 
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Figure 1: The running of the parameters of the soft squared masses a) M? , b) M? 23 , c) 

Mj , and d) M~ when off-diagonal terms are turned off (dashed blue) and on (solid blue) 
at the GUT scale with Aq = 0. Here m 1 / 2 = 520 GeV, mo = 330 GeV, ed = 0.1 and e u = 0.05. 
Also shown is the running when Aq = —ra\j2 and mo = 370 GeV. The dot-dashed (red) curve 
corresponds to the case where off-diagonal terms are turned off at the GUT scale and the dotted 
(red) curve corresponds to the case where they are turned on. 

When Aq ^ 0, the first term of the one loop beta function of ad will now be relevant. 
However since the off-diagonal terms (a/)^ are suppressed with respect to the diagonal 
terms («/)«« as {af)ij/(af)a = e d , for i,j = 2,3 and f — d,e, the difference will be small. 
The running of a^ 2 3 and a e2 3 with A = -m^ is also shown in Fig. [2j 

In order to determine the possible magnification of the flavour violating parameters at 
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Figure 2: As in Fig. [I] for the running of a) a^ 2 3 and b) a e32 . 



the EW scale given the GUT scale initial conditions Eq. (1271) . it will be useful to define 
the parameters, 



0W)« = 7ff=¥=k l (34) 



f,YYJjj 



where (Mj XY )ij, X,Y = L, R are the running soft sfermion parameters. In order to give 
meaning to the flavour violating parameters, we must work in the basis where fermion 
masses are diagonal. The details of this rotation are given in Appendix [B] (see: Eq. ( 1691) ). 
In the definition of the <5's, all soft masses will be assumed to be in the rotated basis. 
In Figure Eh, we show the running of Mj 23 in the rotated basis. Since Mj depends 

primarily on (Y^Yd) 23 ^ e d (if e « ^ e d) when the transformation to the rotated basis is 
performed, we expect that the value of Mj 23 will be reduced due to its alignment with 
the Yukawa matrix Y^. Additionally, for A Q = 0, when off-diagonal terms are turned off, 
there is no contribution from the trilinear terms to the running of soft squared masses, as 
shown in Eq. (1321) . Thus, when going to the basis where Yukawa couplings are diagonal 
(|Y)| 2 = VfYjYfVl*), the contribut ions of the form Y*Yf and YfYj will also be rotated 
away to a large extent. We can see then in Fig. [3^, that this is indeed the case for 
A — 0, ed — 0.1 (blue-dashed line). When A Q ^ 0, we see from Eq. (|32|) that we have 
an extra contribution to the running of Mj from the trilinear terms, which cannot be 
greatly reduced when going to the rotated basis. When off-diagonal terms are turned on, 
for the case of Aq = we expect a greater effect in the running of Mj than in any of 
the other soft squared masses, because (along with Mf) this is the only matrix that is 
sensitive to the parameter at the GUT scale. In Fig. [3h, the (green) long dashed curve 
corresponds to = 0.16 while the (blue) solid curve for = 0.1. Comparing the ratio 
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0.1 2 /0.16 2 = 0.39 we see that in the un-rotated basis (Fig. [Tb) the difference between the 
two initial conditions for M~ will be 39% and is increased to about 50% in the rotated 
basis because the diagonalization matrix also contains the parameter e d . 

The sensitivity to e u in the running of Mj is very small as a result of the difference 
between e 2 , and e 2 . The sensitivity to c 22 in the running of the soft squared masses comes 
from the terms of the form YfYj and aja^, however as we will describe in the next section, 
the elements of Yy . . are adjusted to reproduce the correct fermion masses at the EW scale 
for a given choice of c 22 and e^. For c 22 = 0.7 and = 0.1, one is able to reproduce the 
correct value of m s (EW). Similarly, for c 22 = 0.5 and = 0.12 we are also able to obtain 
the correct value for m,(EW). Indeed, the product c| 2 e^ at the GUT scale is very similar. 
Therefore, the quantities YfYj are more or less fixed by this requirement for different 
choices of c 22 and e 2 d and result the running of Mj will not be very sensitive to these 
changes. 

In the remaining panels of Fig. [3j we show the flavour violating parameters (5 qt xx)23, 
(8cI,xy)23, and (<5 e ,xx)23 for our test point (m 1 / 2 , m , A ) = (520,330,0) GeV and 
(erf, e«) = (0.1, 0.05). The curves for these cases are all shown as solid (blue). Off-diagonal 
terms at the GUT scale are turned on for all curves in panels b) - f). For comparison, we 
show the case with non-zero A , ie., (my2, trio, Aq) = (520,370,-520) GeV. This case 
is shown by the dashed (blue) curves. In addition for both A = and A ^ 0, we show 
the sensitivity of the running to e^. The (red) dotted and dash-dotted curves correspond 
to these cases respectively, where we have chosen = 0.16. 

6.2 Determination of coefficients of Yukawa matrices 

In this subsection, we will determine the values and sensitivities of the quantities 

C 23> C 33> C 22? C 23' C 33> C 22 an d c 33 (35) 

to the variation of and e u for a given value c 22 . For now, we continue to use our 
example based on the point (j7ii/ 2 , m , A ) = (520,330,0) GeV, tan/3 = 50. As a result 
of the small differences in the running of the off-diagonal terms of soft parameters when 
we choose different initial conditions at the GUT scale, the Yukawa parameters are not 
really sensitive to this difference. In this section, we will only consider cases where the 
off-diagonal elements are turned on at the GUT scale. These parameters are, however, 
sensitive to variations in e d , e u and c 22 . 

In Fig. HI we plot the relevant qj- versus the appropriate e/ for c 22 = 0.7. There is 
little sensitivity to td or e u for the 33 components of the Yukawa matrices which are well 
fixed by the low energy fermion masses. In contrast, the 22 and 23 components can be 
sensitive to and e u even within their limited range of variation. 

In Fig. |5l we show the dependence of two of the Yukawa parameters, c 23 and cf 3 , to 
our choice of the supersymmetric point as a function of e^. In addition to our original 
test point (shown here as the dotted curves), we consider (m^, m , A ) = (160,630,0) 
GeV shown by the dashed curves, and (m^, m , A ) = (1240,910,-1240) GeV shown 
by the solid curves. In all three cases, we have chosen tan (3 = 50. In each case, the value 
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Figure 3: In panel a), the running of Mj 23 in the rotated basis as in Fig. [TJ Also shown 
by the (green) long-dashed curve is the case = 0.16 with Aq = and off-diagonal terms 
turned on at the GUT scale. In the remaining panels, the running of squark flavour violating 
parameters (5 q ,xx)23 (b-d), (5 djX Y)23 (e), and (5 et xx)23 (f) for (m 1/2 ,m , A ) = (520,330,0) 
GeV, with tan/3 = 50, = 0.1 and e u = 0.05 shown by the solid (blue) curves. In panels b) 
-f), off-diagonal terms are on at the GUT scale. Also shown by the dashed (blue) curves are 
the results for Aq = -mj/j with wiq = 370 GeV. The (red) dotted and dash-dotted curves show 
results for = 0.16 with Aq = and Aq = —520 GeV respectively. 



18 




Figure 4: Yukawa parameters, c£ of the d, e and u sectors for the supersymmetric point 
K/2,mo, A ) = (520, 330, 0) GeV, tan/3 = 50 and c% 2 = 0.7. In panel (a) we show cf 3 for i = 2, 3 
as labeled as a function of for fixed e u = 0.05. In panel (b) we show cfj for ij = 22, 23, 33 as 
labeled as a function of for fixed e u = 0.05. Finally in panel (c) we show for i = 2, 3 as 
labeled as a function of e u for fixed = 0.1. 



of mo was chosen so as to obtain a relic density compatible with WMAP. In the latter 
case, we chose A = -m 1( / 2 . The dependence of the parameters of Y e on the change of 
(mi/2, mo, Aq) is similar though less pronounced and the dependence of the parameters of 
Y u is very small. 

Given the parameterization of the Yukawa matrix in Eq. (|22|) . the choice of c 22 is 
directly correlated to the allowed value of given the value of m s /mf,(MGUT): 

m s (M \ (^22 (j^)j2 + (j^)g3 C 22 2 4 

^ (Mgut) ~(T^ 2(5?)^— -S"" 0( °' (36> 

which must be in the range ~ (0.01, 0.03) at Mqut- Hence for different choices of c! 
assumed to be 0(1), we have the following allowed ranges of ej. 



d 

22' 



4, = 0.5 -> e d ~ [0.10,0.12] 
4j = 0.7 -> e d ~ [0.085,0.10] 

4, = 1.0 -> e d ~ 0.08. (37) 

As a consequence the values of (Yd)ij and (Y^.)y are different for each case. In Table (TJ we 
show the allowed ranges of the Yukawa parameters according to the allowed range for 
in each case. 

To further compare how the determination of the appropriate Yukawa parameters at 
the GUT scale depends on the chosen supersymmetric parameters we show in Fig. [6j the 
values of c£ as a function of the the gaugino mass, m 1 / 2 - Here A = and tan/3 = 50. 
For each value of m^, m is adjusted to insure the correct relic density of neutralinos. 
As one can see the dependence on the supersymmetric model is rather minimal. 

Finally, we would like to emphasize the importance of taking into account the su- 
persymmetric corrections to the fermion masses for the proper determination of Yukawa 
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tanP = 50, e u = 0.05 



0.08 0.09 0.1 0.11 0.12 0.08 0.09 0.1 0.11 0.12 

Figure 5: As in Fig. HI we plot (a) and (b) C33 for three different choices of (my 2 , mo, A)). 
The bottom curves (dotted) correspond to (mj/2, mo, Aq) = (520,330,0) GeV, the lines in 
the middle (dashed) to (mi/2, mo, Aq) = (160,630,0) GeV and the lines at the top (solid) to 
(m 1/2 ,m ,A)) =(1240,910,-1240) GeV. 





Yukawa parameters 






c 22 


c 23 c 33 c 22 


C 23 


C 33 


0.5 


[0.6,0.9] 0.3 [1.2,1.7] 


[1.3,1.7] 


0.5 


0.7 


[0.9,1.3] 0.3 [1.7,2.4] 


[1.9,2.4] 


0.5 


1.0 


1.6 0.3 3.0 


3.1 


0.5 



Table 1: Values of the Yukawa parameters for different choices of C22 for the supersymmetric 
point (m 1/2 ,m , A ) = (520,330,0) GeV. 

parameters at the GUT scale. From [39], we find that for large tan/3, the supersymmetric 
corrections to the mass of the bottom can be as large as 30%, as it is for the case of 
tan/3 = 50. We can express the relation between the masses of the fermions at a certain 
scale Q, containing the supersymmetric corrections, and the masses of fermions at the 
same scale as given by the SM through: 

m s k (Q) = m s k M (Q)gl(Q), k = t,b,c,s. (38) 

Here mf. M (Q) are the SM values at Q, including the corrections from gluons. A good 
analytical approximation to the factors g k {Q) is given in [39] . For example, for the 
top quark, this can be expressed as gf(Q) = 1 + Am t /m t , where the most important 
contributions to Am t /m t comes from top squark/gluino loops. In Fig. we plot the 
numerical values of g k (Q) for k = t,b,s, as a function of m\i% for tan/3 = 50 and for 
the choice of Yukawa parameters (c^q) = (0.7,0.1). These corrections include the 
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Figure 6: The coefficients of the Yukawa matrices plotted versus m X j 2 for = 0.1, e u = 0.05 
and C22 = 0.7. In panel (a) we show cf 3 for i = 2, 3 as labeled. In panel (b) we show cfj for 
ij = 22, 23, 33 as labeled. Finally in panel (c) we show for i = 2, 3 as labeled. 

corrections at one loop from squark/gluino loops, charginos and neutralinos. We can see 
that, in fact, for m& these corrections are about 30%, while for m s , they can be up to 
14%. For m t , the corrections are within 11%. The corrections for m c , though included 
are not shown as they are very small. 

The supersymmetric corrections change the profile of Yukawa matrices at the GUT 
scale. A SM analysis [2] determined the values for Yd (for symmetric matrices) as follows: 

= 0.15 2 = 0.0225, = 1.3(0.15 2 ) = 0.0293. (39) 

(Yd) 33 [Yd)33 

If we take the values of the Yukawa parameters for (c^, q) = (0.7, 0.1) from Tabled], we 
have 

^hl = 0.025, = 0.045. (40) 

(Yd) 33 (Yd)33 

Since c% 2 is an input to our analysis, we are basically adjusting the value of (1^)22/(^)33, 
at the GUT scale, thus there is no surprise in the similar values for this ratio in both 
cases. What is interesting is the change in the off-diagonal parameter (1^)23/(^)33, since 
this off-diagonal term is crucial when constructing a specific flavour symmetry. 

It is important to emphasize the difference between the structure of the Yukawa ma- 
trices at the GUT and EW scales. Whilst at the GUT scale they can be written in terms 
of the expansion parameter e 2 d and the prefactors c{j, at the EW scale this is not longer 
the case. However if the prefactors qj- are indeed to be explained by a specific FS then 
the result at the EW can be expressed as a series expansion in terms of e/ and the same 
coefficients. For example, we can express the Yukawa matrices normalized to the elements 
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Figure 7: MSSM corrections at one loop to the masses of m t (solid blue line), m& (dotted green 
line) and m s (dashed red line). 



(17)33 at the GUT scale. The numerical value using c% 2 = 0.7, = 0.1 and e u = 0.05 is 
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where the effect of running to the EW scale is given to the right of the arrows. 
6.3 Effects on b -> s 7 

We have seen that a good choice of Yukawa parameters was c% 2 = 0.7 and = 0.1, 
e u = 0.05, thus we present the analysis on the sensitivity of b — > s 7 on centered on 
those parameter values. As we explain in Appendix El due to the hierarchy of <C e|, we 
expect that the main departures from MFV come from the gluino contributions, entering 
through the magnetic Wilson coefficients due to the diagonalisation of the soft squared 
matrices in the d sector, Eq. fl27j) . In Fig. [8] we show our calculation of b — > s 7 as a 
function of when the running of the parameters of the second family have not been 
considered (dashed blue curve). For comparison, we also show the same result when the 
running of both the 2nd and 3rd generations are included for Aq = (solid red curve) and 
when A = —m 1 / 2 (dotted black curve). As before, for each value of m 1( / 2 , m is chosen so 
as to obtain the WMAP value for the relic density of neutralinos. Note that the inclusion 
of the 2nd generation running relaxes slightly the bounds obtainable from a comparison 
of the computed value of b — > s 7 with the experimental measurement at relatively low 
Tii/2 (~ 400 GeV). As long as we fix the relation between C22 and through the fermion 
masses, there will be little sensitivity of b — > s 7 to e^. 
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Figure 8: The branching ratio for b — > s 7 for Aq = (solid red) and Aq = —mu2 (dotted 
black) when e<j = 0.1, = 0.05 and C22 = 0.7. When only the running of the 3rd generation is 
included, the analogous result is shown by the dashed (blue) curve. 

The contribution of the charged Higgs to the Wilson coefficients, depends crucially on 
x t H = ml/m 2 H - where m#- is the mass of the charged Higgs. The determination of fun- 
is quite sensitive to changes in the running of the supersymmetric spectrum [ID] and thus 
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any change in the running of these parameters will affect the final value of m#- : 

m 2 _ 2 , M 2 , A 2 A 2 _ -302 m tV 2 /( M Ii,« 2 ) r42 N 

where is the pseudo-scalar Higgs mass and /(Mj fi ) is a function of the lightest sups: 
Ui for z = 1,2. The running of /i is quite sensitive to the changes in the running of the 
susy spectra and as a result m 2 H - is also sensitive. Thus, when we include the running 
of the second family, the added small changes to the supersymmetric spectrum alters the 
final value of tuh- ■ 

K 2 

We have tested the effect of the gluino contribution by computing the ratio, 5,~ " — - J ~ d ~ 



g k 2 ' 

7offg 

where these quantities have been defined in Eq. (1551) and Eq. (1631) in Appendix |A] The 
ratio is very close to 1 for mi/2 > 300 GeV. At lower gaugino masses, the contribution 
from gluinos increases but the ratio only decreases by about 1% at m\i% = 160 GeV. 
When mi/2 is below 300 GeV, the stops become quite light and a NLO calculation is not 
possible using the approximation of [TTj . At lower mi/2, a LO calculation has to be used 
also for charginos, charged Higgs and neutralinos. In fact when < 280 GeV this is 
the case. Thus all the results for b — > s 7 constrained to < 280 GeV are calculated 
at LO. 



6.4 Constraints on the parameter 

As we have pointed out in Section HI once a family symmetry is introduced, additional con- 
tributions to the trilinear couplings are given by the derivatives of the Yukawa couplings 
with respect to the flavon fields, ~ Y/jF a d a , and then in a realistic flavour symmetry 
{a>f)ij must be of the form of Eq. (TT6T) with kL 7^ 0. Given the limitation on the number 
of parameters that we can use to scan the possibilities of the general form of Eq. (f2"Tl) . 
with coefficients x{, 7^ 1 and z{j 7^ c{j, we choose to analyze the matrices of the form 
given in Eq. (1281 . and hence effectively allowing e$f > e/. To remain consistent with the 
approximations made earlier, we maintain the restriction es u < 

To this end, we first check once more the running of the off-diagonal squared masses 
and trilinear couplings with respect to the parameters e S j > e/. For the case of e S f = ej 
we have seen in the discussion of section [6TT1 (Fig. [[]) that (M|) 23 and (M|) 2 3 are more 
sensitive than the quark sector, other than (M|) 2 3- To be definite, we make the following 
choices for esf- 

/■ e Sd = e Se = 0.45, e Su = 0.06, 

II- £sd = 0.45, e Se = 0.225, e Su = 0.06, (43) 

for both cases we keep c\ 2 = 0.7 and e d = 0.1, e u = 0.05. In Figure El we plot then the 
running of M| 23 , Mj 23 , Mj 23 and M| 23 for the two cases of Eq. fT43l . These figures can 
be compared directly with those in Fig. [lj The only difference being the increased values 
of the expansion parameters which are now distinct from those in the Yukawa sector. 

The running of Mq 23 is significantly affected by the increase in e sd as one can see from 
the comparison between Figs. [T^, and [9h. This is largely due to the increased importance 
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Figure 9: As in Fig. [TJfor case I with e$d = tse = 0.45, es u = 0.06 (red dashed lines), and for 
case II with egd = 0.45, es e = 0.225 and es u = 0.06. (blue solid lines). Here mo = 370 GeV and 
A = -m 1/2 . 



of the a d a,d term in the beta function which is now competing with the Yukawa terms in 
Eq. [2HJ Not surprisingly the running of Mj 23 is quite close to that of Figure [Qd because 
we have kept es u = 0.06 quite close to es u = 0.05. In order to affect the branching ratio 
for b —>■ s 7, we are most interested in the effect on the d sector. The initial value of 
Mj 23 is determined by esd (and mo) and is increased relative to the case studied in Fig. 
[Tb by the the ratio esd/^d = 0.45 2 /0.1 2 = 20.25. As a result, we can see in Figure the 
enhanced running of Mj . In panel d), we see the differences between cases I and II, as 
the running of M| 23 is clearly sensitive to es e - In Fig. [10J we show the corresponding 
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IjhxYUl for (m 1/2 ,m ,A ) = (520,370,-520) GeV 







I 


II 


\( S L,ll)w\ 


0.0012 


0.15 


0.014 


(4, M )23 


0.0069 


0.50 


0.05 


(4, M )23,32| 


6.0 x 10~ 6 


8.6 x 10~ 3 


1.5 x 10~ 3 


(4 M )23,32| 


4.4 x 10~ 6 


7.2 x 10~ 3 


1.1 x 10~ 3 


\(^Q,ll)^\ 


0.0069 


0.009 


0.009 


\(fid RR )23\ 


9 x 10" 4 


0.021 


0.021 


(Vla)23 


4.3 x 10" 6 


8.0 x 10~ 4 


8.0 x 10- 4 



Table 2: Weak scale values for the |(<5/,xy)23| (denoted in the 1st column) for (mi/2,mo,Ao) = 
520, 370, —520) GeV. In the 2nd column, we assume esd = = 0.1, esu = e« = 0.05. In column 
three, we take esd = 0.45, es u = 0.06 and = 0.1, e u = 0.05 while in column four, we lower es e 
to 0.225. 



effect on the trilinear couplings which can be compared directly with Fig. [2j 
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Figure 10: As in Fig. [9] for the running of a) 0^23 and b) a e32 . 

In Fig. [TTJ we show the effect of the increased es on the parameters defined 
in Eq. (jMl) . The increased effect on the running of these parameters is also largely due to 
the increase in the initial values. For example, for (^,7^)23 the initial value is increased 
from 9 x 10~ 3 to 0.021, when = egd changes from 0.1 to e$d = 0.45. The comparison to 
the increase for the other (5/ i xy)23 is given in Table [2J 

If one would like to introduce right handed neutrinos for these cases, then the present 
analysis will be valid assuming the right handed neutrinos decouple at a higher scale than 
^gut, or M UR > Mqut- For M VR < Mqut the quark sector will not be largely affected 
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Figure 11: As in Fig. [9]for the running of (Sf t xy)%j with /, XY = a) e, RR b) d, RR, c) Q, LL, 
d) u,RR, e) d,LR, f) e,LR. 

but the lepton sector will be. We will have then important changes to the running of 
(M|) 23 and (Mf ) 23 . Since the goal in this paper is to study the constraints mainly in the 
quark sector, the introduction of M UR < M GU t is beyond the scope of the present work. 
For the present case we consider constraints from leptonic flavour violating process such 
as t — > ytry. For example in the CMSSM, the bounds in terms of the flavour violating 
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parameters Eq. (|34"|) . for tan/3 = 50, are [4"T] 



|(5L,LL)23| < 0.1, \(5 e>RR )23\ < 0.1, 1(^)23,321 < 0.02, 



(44) 



using the most stringent current bounds. For the case es/ 
($i,xx)23 are well below these bounds: 



€f the values obtained for 



\{^L,LL)2?\ 
|(^e,L_R)23| 



0.0012, 
6.0 x 10" 



|(5 e> _RL)23| 



0.0069, 
= 4.4 x 10" 6 . 



(45) 



For case I of Eq. fT43|) the bounds of Eq. (JHJ) are already saturated while for case II, the 
parameters are still within the experimental limits. 

As can be expected, the sensitivity of b — > s 7 on our expansion parameters for the 
cases €ds > e<f is greater than that discussed in the previous section with ejs = The 
value for the branching ratio for b — > s 7 as a function of m 1( / 2 is shown in Fig. [121 
The lower two curves correspond to the single family case (where all mixings are ignored) 
and our nominal choice of e$d = Q and are taken directly from Fig. [HJ The upper two 
curves show the effect of increasing esd and es e using values from case I (upper solid red 
curve) and case II (green dot-dashed curve). Here we see more clearly that in regions 
with relatively low m^, the effect of mixings relaxes the constraint imposed by b —>■ s 7. 
In these cases, the gluino contributions can be up to 10% for mi/ 2 around 200 GeV but 
when mi/2 increases these contributions reduce considerably. Here the contribution from 
gluinos does not vary too much from cases I and and II of Eq. fT43|) , as the running in the 
d sector is not greatly affected. 
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Figure 12: As in Fig. [H The curves from bottom to top are as follows: the single family limit 
is shown by the blue dashed curve; the choice of = e^^ = 0.1 is shown by the black dotted 
curve; increased mixing with esd = 0.45, es e = 0.225 is shown by the green dot-dashed curve; 
and e$d = ^Se = 0.45 is given by the red solid upper line. 
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To better appreciate the effect of the mixing in the squark sector, we zoom in on 
the region between muz = 300 and 800 GeV as shown in panel b) of Fig. [TJJ For 
tan /3 = 50, this region is preferred when other observables such as (g — 2) M are included 
[33l W2\ . When theoretical and calculational errors are included, the 95% exclusion value 
for BR(6 s 7) occurs at about 1.9 xl0~ 4 , as discussed in the previous section and 
corresponds to a lower limit of 400 GeV on m^j in the single family case most commonly 
treated in the CMSSM. This limit is substantially relaxed when two-family mixing is 
included. For es e = £sd = 0.45 as shown by the upper curve in the figure, we see that the 
limit is relaxed to > 335 GeV. One should keep in mind that the exact numerical 
limit found here and in the CMSSM is subject to change when NNLO SM corrections are 
included. On the other hand, we expect the relative change in the limit when two-family 
mixing is included to be a robust result. 

Finally, we show in Fig. [T3] the behaviour of BR(6 — + s 7) on esd f° r several choices 
of mi/2- As one can see, for low mj/2, the branching ratio is lower than the experimental 
value when esd is small. When esd > 0.4, the branching ratio increases and the constraint 
on mi/2 is relaxed. 
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Figure 13: Variation of BR(b -t s 7) with respect to esd- From top to bottom the lines 
correspond to (m l/2 ,m ,A ) = (600,410,-600) GeV, (500,360,-500) GeV, (400,340,-400) 
GeV and (300, 380, -300) GeV. 



6.5 Sensitivity on the determination of the MSSM spectra 

In this last section, we would like to consider the effects of the our flavour violating 
parameters on the resulting supersymmetric spectrum. Since the parameters e u , e^, es u , 
£se, £sd are small (< 1) and their dependence on the Yukawa, trilinear and soft squared 
matrices is quadratic, Eqs. (I27H28|) . we do not expect a big effect on the spectra, except 
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for sensitive parameters such as 



fi, m Hl2 , rriA, 



(46) 



which are determined from the electroweak symmetry breaking conditions. In Fig. [Tj 
we show the dependence of the six masses of the sdown sector and of the charged slepton 
sector with respect to e Sd for the case 



;£sd, 



(47) 



with esu = 0.06 which includes case II of Eq. ()4"3l) for e Sd = 0.45. Here we have fixed 
m l / 2 = —A = 520 GeV, m = 370 GeV, and tan/3 = 50. As one can see, only when 



tan B = SO , m 1/2 = 520 GeV, A = -520 GeV tal > P = 50 , m 1/2 = 520 GeV, A = -520 GeV 
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Figure 14: Masses for the sdowns (a), and charged sleptons (b) sectors as a function of e$d for 
the case (mm, mo, -Ao) = (520,370,-520) GeV, es e = £d = 0.1, and e u = 0.05. In panel c, 
the soft Higgs masses, \i and txia are shown as labeled. 

e sd ?t 0.35, does the splitting in the masses of the heaviest families (1st two generations) 
become relevant. The third generation (mostly sbottom) masses are relatively insensitive 
to tsd and are close to their CMSSM (with no flavour violation) values (Mj i? Mj 2 ) = 
(956,1017) GeV as shown by the lower two (red) curves in Fig. [T4"h . where 1,2,. ..5, 6 
are ordered lightest to heaviest. Analogously for the 3rd generation sleptons we have 
(Me 15 M e - 4 ) = (246,477) GeV, in the CMSSM as shown by the two (red) curves in panel 
b. Note that the heavier stau is heavier than the mostly right-handed selectron and smu. 
With no flavour violation, the masses of the 1st two generations are degenerate (due to 
the smallness of the associated fermion masses) in the CMSSM. Here we see the degree 
of splitting induced when esd ^ 0. We do not show the dependence of the of the sups 
with respect to esd or es M because their change (with respect the running of just the 
third family) is less than the 0.4%. We note, however, a more important difference in 
the determination of the running Higgs sector parameters which is displayed in panel c 
of Fig. H3J When e Sd = 0, we have (ji,m Hl ,m Ha ,m A ) = (568,-403,-708,516) GeV at 
the electroweak scale, where m# 4 are the two soft Higgs masses, and the minus refers 
to the sign of the squared mass. For esd = 0.45 these become (/i, m^, m# 2 , uia) = 
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(545,-465,-688,452) GeV. Since the most important supersymmetric contribution to 
BR(B S — > scales as /§ /m\ and in the analysis that we have done the inclusion of 

the second family tends to lower m^, we expect that in general BR(£> S — > increases 
from the case egd = to e$d = 0.45, but remains within the experimental limit. 

7 Summary 

Flavour symmetries (FS) are often constructed using supersymmetry. Departures from 
MFV will affect the supersymmetric spectra and constraints at the electroweak scale 
can be employed to determine the profile of the Yukawa matrices at the GUT scale. 
Flavour violating processes then can be used to further constrain the shape of such Yukawa 
matrices and the parameters associated with the FS. 

Clearly, a completely general analysis is not practical for computational purposes, as 
it involves many new supersymmetric parameters in addition to the Yukawa matrices. To 
constrain the parameters of the latter we can rely on the determination of the Yukawa 
matrices at Mew from flavour violating parameters in the SM, assuming a hierarchical 
symmetric form, Eq. ([6]) (although in this context a non-symmetric form is equally plau- 
sible, it requires more parameters) such that the expansion parameter in the up sector is 
smaller than in the down sector, motivated by the different hierarchies between up and 
down quarks respectively, e u ~ e 2 d . For the latter case we can simplify the usual rela- 
tions obtained in FS between the expansion parameters of the Yukawa matrices and the 
off-diagonal parameters of the trilinear couplings and soft squared mass matrices at the 
GUT scale. 

In our analysis, we used an iterative procedure as is common in the analysis of the 
CMSSM [28| [29] . with the inclusion of three more parameters describing the Yukawa 
matrices (e M , and c% 2 ), at the GUT scale. Other necessary inputs are the values of 
m s (2GeV), m c (2GeV), m M (M EW ), and the values of the CKM matrix relevant to the two 
heaviest families. These are used to determine seven more parameters in the Yukawa 
sector at the GUT scale, as in Eq. f[2"2"j) . 

The important results of this analysis in the fermion sector are as follows. First, us- 
ing the supersymmetric corrections to the fermion masses of the two heaviest families, 
at the EW scale, we can accurately determine the profile of the Yukawa matrices at the 
GUT scale. Particularly important is the ratio (Y rf )23(M GUT )/(Y rf ) 2 2(M GUT ). Whereas in 
an analysis where the running of supersymmetric parameters is neglected, this ratio is 
found to be 1.3, when supersymmetric corrections are taken into account, this ratio is 
enhanced to 1.8, cf., Eq. fT40l . Second, the value of the expansion parameter, e^, is effec- 
tively constrained by the appropriate determination of m s (MEw), with its corresponding 
supersymmetric corrections. Of course the exact value depends on the values for the 
coefficients used in the expression of the Yukawa matrices, for example (Y d ) 2 2 = c^e^. 
Assuming that the coefficients cfj should be of order one, we find for example that for 

= 0.7, td G [0.085,0.10], otherwise m s is predicted to be either too small or too big, 
respectively. This determination is important when constructing FS. The exact values of 
(1/);,,(M G ut) are used to justify different choices of family groups (e.g. SU(3), SO(3)). 
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We have only considered the running of the two heaviest families and hence the structure 
of the lightest family cannot be determined at the GUT scale by this analysis. However, 
the inclusion of the lightest family will not have a strong impact on the results obtained 
here. 

Our results are also relevant for the supersymmetric sector. When off-diagonal terms 
are introduced in the sfermion mass matrices or trilinear terms at the GUT scale, de- 
partures from MFV are induced at the EW scale. In the most constrained models, ie., 
when no new parameters are introduced as in Eq. ( 1271) . this departure is small and the 
resulting sensitivity of BR(b — > s 7) is not significant. In contrast, when the expansion 
parameters, esf , for the supersymmetric sector differ from that of the fermion sector, e/, 
we may obtain more sizable departures from MFV. We have compared the running of the 
off-diagonal trilinear couplings and soft squared mass matrices for these two possibilities. 

Our calculations were performed using the values of the supersymmetric parameters 
as obtained when running the full MSSM, with the simplified structure for the trilinear 
couplings and soft squared masses at two loops as described above. Hence we do not rely 
on estimating the value of the flavour violating parameters (Sf XY ), Eq. (134|) . at the EW 
scale. Indeed, these are calculated directly in terms of the parameters introduced at the 
GUT scale as an extension of the CMSSM. These parameters were introduced in Section 
16. II to emphasize their running from Mqut to Mew- Current limits on the charged slepton 
sector could be used to further constrain the choice of the expansion parameters e/ and 
e sf- 

We have seen that the running of (Sf XY ) 2 3 is quite limited when initial conditions at 
the GUT scale are constrained. This is true for both the case of MFV at the GUT scale as 
well as the case when the expansion parameters of the supersymmetric sector differ from 
those of the fermion sector, eg/ 7^ e/, as we saw from Table [2j In column 2 of that table, 
we have assumed MFV at the GUT scale and the departures from MFV at the weak scale 
are small as seen by the small values of the ^-parameters. In this case, the bounds on the 
parameters (5( e ,z,) xr )23 are well below the bounds imposed by BR(r — > fij) in Eq. (1441) . 
However, if we take for example e$d = £se — 0.45, we can see that the values obtained for 
($(e,L) XY )23 in column 3 already saturate the bounds imposed by BR(r — > jj,^). In the last 
column, we showed the 5's when esd = 0.45 and es e = 0.225 where the bounds imposed 
by BR(r — > fi-y) are satisfied. 

We have linked the down-type quarks to the charged leptons by generating the off- 
diagonal terms of the Yukawa matrices of both with Vckm at EW scale as in Eq. f[2"3"j) . This 
is motivated by GUT theories where often charged leptons and down-type quarks acquire 
Yukawa matrices of the same form. Extending the relation to the supersymmetric sector, 
by taking the expansion parameters esd and es e to be of the same order, we naturally 
find that the parameters \(8{Q,d)xY)2z\ and |(£(L, e )xy)23| are correlated at EW scale. From 
Table El we see that for e Sd = e Se we have in fact \{5 {Q4 )xx)23\ = O (\(5( Lte ) X x)zs\)- 
Applying the bounds on |(5(£, je )xx)23| imposed by BR(r — >■ /ij), we have taken es e = 
0.5esrf = 0.225 and this effectively decreases |(5(£, e )xx)23| by an order of magnitude with 
respect to K^q^xx^I- 

Finally, we have calculated BR(b — > s 7) at NLO for the charged Higgs, chargino and 
neutralino contributions and at LO for the gluino contributions (since a NLO computation 
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is not yet fully available for this case). BLO corrections for the gluino contributions were 
considered in [TS], but given the values of (<5/ xy )23 that we have derived, these corrections 
will have little impact on the value BR(b — > s 7) that we obtain. 

Given the SM NNLO computation and the continuously improved experimental value, 
it is important to determine BR(b — > s 7) as precise as possible in a given supersymmetric 
model. We have seen that the structures for Yukawa matrices, trilinear terms and soft- 
matrices of the form of Eqs. (127H28I) yield a positive contribution to the calculation of 
BR(b — > s 7) for fi > bringing the branching ratio closer to its experimental value than 
what would have been obtained in the absence of any FV. In fact when esf 7^ £/, such that 
£s<i > Cd > e s« > e„, it is possible to relax the lower limit of ~ 400 GeV on mi/2 found in 
the single family case most commonly treated in the CMSSM. For es e = tsd — 0.45 this 
limit is relaxed to myz > 335 GeV, as we have discussed in Section |6"41 



A b — > S7 beyond minimal flavour violation 

At LO the renormalization evolution of the gluino contribution to the effective Hamil- 
tonian, H e s, from Mew down to //& is different than the contribution from W~ , H , 
charginos and neutralinos [Hj. Hence, TC e g can be written as: 

H eS = H™-' H -'* + Hl s . (48) 

"Mixed" diagrams containing the W boson, gluinos and squarks give rise to a s corrections 
to the Wilson coefficients in 7i w < H ' x at the matching scale, and they are taken into 
account at the NLO level. We include the contributions of gluinos at LO and those from 
H~, charginos and neutralinos at NLO. 

We recall that the starting point in the calculation of inclusive B decay rates is the 
low-energy effective Hamiltonian [20] 



AC 

H ~' X = -^TOE^O^) . (49) 
In the massless strange quark limit, the operators relevant to our discussion are 

lbir z 

Os = g -^s L a^t a b R . (50) 

To an excellent approximation, the contributions of other operators can be neglected. 
The renormalization scale /z& in ( 14*91) is conveniently chosen of order m^, so that all large 
logarithms reside in the Wilson coefficient functions. However when the mass of the 
strange quark is taken into account we need to consider the operators 

0' 7 = ^s R a» u F^b L , 
167r 

O's = g -^s R a, u G^t a b L . (51) 
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For these operators the corresponding Wilson coefficients C H oc m s mb/m 2 tan 2 (3 and 



C x are quite small. For TCl s we have 

Kn = Y, c uM°<aM + E^(^)o«(^) + EE c 2!ff(^)Q®(^). ( 52 ) 

i i i <1 

where q runs over the light quarks q = u,d,c,s,b. The operators can be of 

two types, one for Oj bjj{.Hb) and the other for Oj g,g(fib)i j = 7,8. The former can 
be found from Eq. ( 1501) with the substitution of (1/167T 2 ) with g 2 and the latter with 
the substitution of (m;,/167r 2 ) with g 2 . The primed operators in Eq. (1521) come from the 
interchange of L <-» i£. The operators Of§(/^) are suppressed with respect to the magnetic 
operators by the mass of the corresponding quark or by an extra power of g s . We neglect 
them here since their contribution is quite tiny with respect of the operators in the first 
two terms of Eq. (j52|) . 

In the SM, the complete theoretical prediction for the B —>■ decay rate at NLO was 
presented for the first time by Chetyrkin et al. [21] . It depends on a parameter S defined 
by the condition that the photon energy be above a threshold given by > (1 — 5)E™ ax , 
where = mb/2 is the maximum photon energy attainable in the parton model. Thus, 
5=1- 2£™ ax /m fe , and we use £™ ax = 1.6 GeV, as it is the value at which the HFAG 
collaboration has presented its latest results [HI EH] • 

The prediction for the B —>■ Xgj branching ratio is usually obtained by normalizing 
the result for the corresponding decay rate to that for the semileptonic decay rate, thereby 
eliminating the strong dependence on the 6-quark mass: 

BR(B -> X s 7)^ 7>(1 _ 5)K nax = BR(B -> X c eu) cxp — — — — - — ^— . (53) 



Following closely the formalism of [221 123] , we write 

v t * s v tb 



BR(B -> X s1 )\e i>Eo = BR(B -> X c ev\ 



exp 



V, 



cb 



2 6a. 



cm 



T(B -> X c eP) 



[P(E ) + iV(E )] - /T ,, (54) 



7T ' r(r u /r c ) 



where perturbative, P(E ), and non-perturbative, N(E ), contributions are treated sepa- 
rately, and Eq = (1— 5)E™ &X . The function N(E ) is a correction which does not depend on 
the renormalization of the running mass from /xew to m^, it contains power corrections 
that arise from the interference between the current-current operator and the magnetic 
dipole operator, hence called the "non-perturbative" contribution. The precise definition 
is given in Eq. 3.10,ofref. [22]. In the SM this is [22], N(E = 1.6 GeV) = 0.0036±0.0006. 
The ratio r(r M /T c ) is a semileptonic phase space factor which is independent of the cal- 
culational problem of convergence of the perturbation series in b — > Xgj. This has been 
calculated accurately [22] up to NNLO in the SM. 

Since BR(B — > X c e^) is an experimental input, only the function P{Eq) has contribu- 
tions from processes beyond the SM. Using the approach of Gambino et. al. [22], rather 
than writing the terms P(E ) + N(E ) as functions of the standard Wilson coefficients, 
Eq. (jUP , at the decay scale, fib, these functions can be written in terms of functions 
involving the renormalization effects from /i (m t or /xew ) to \i\> and the Wilson coef- 
ficients at fiQ. The contributions from the charm quark do not get re-normalized from 
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Mew to m b , while the contributions from t do. Hence instead of the familiar contribution 
of \Cf (/i E w)| 2 to BR(B -> X s j), we have 

|Cf(/i fe )| 2 - K 2 7cS = P(E ) + N(E ) 

P(E Q ) = \X c + X t + e EW \ 2 + B(E ), (55) 

where X c and X t represent the c and t contributions respectively, €ew is an electroweak 
correction and B(E ) is the bremsstrahlung correction from b — ► 57^ and 6 — ► S7gg with 
q = u, d, s [23J. 

Analogous to the usual Wilson coefficients, we can decompose 



X/M = Xj 0) (/i fe ) + Zj^xVfa), f = c,t (56) 
The effects of physics beyond the SM can be incorporated [23] as: 

Cf> T0T (/, EW ) = Gf )SM (/iE W ) + cf "^W), < = 7, 8, (57) 

where C|°^ SM (//ew) — C ( ' ' ,efr (A t Ew)i, i = 7,8 correspond to the coefficients of Eq. (12) 
of [21], i.e. the scheme independent Wilson coefficients, the coefficients Cf^ BSM (/^Ew) 
of course refer to processes beyond the SM, in this case the MSSM. From here on, we 
will refer to the following coefficients C 7 °\h-ew) = C 7 °^ TOT (/jew) + §§ an d Cg (#ew) — 

^ 0)TOT (/i E w) + |. For the MSSM, the general form for Cf )MSSM (/i EW ) is known, al- 
though the contribution from gluinos does not follow the same pattern of renormalization 
from Mew to mj and it should be implemented in a different fashion. The contribution 
from charginos, neutralinos and charged Higgs bosons can be renormalized as in the SM, 
although the complete NLO corrections to (7^| s 8 M (/i{,) are not known. 
Thus, if we have the following hierarchy of masses 

/isusv ~ 0(rrig, m$, m^, m t ~ 2 ) > /i E w ~ 0(M EW , rn t , m H ±) > fi b , (58) 

we expect that the following assumptions represent a first good approximation for taking 
into account the contribution from charginos, neutralinos and charged Higgs bosons: 

Cf = CP SM + C^ H± , i = 7,8. (59) 

We note that C 7 (n b ) must be added on directly to xj: \u, b ) because gluinos undergo a 

different renormalization from Mew to \i b . The expressions that we use loiCfl x ' H and 

Cf ] are given in Eqs. A16-A20 of gU- At NLO for (7* H± , we follow the approach of [Hj. 
Thus at LO, Eq. (155)) would receive the following contributions from gluinos at //&: 



P(E ) g = xi 0) (ji b )+X c (ji b )+X t (ji b ) + e EW ~+ Xf (v b )- + B(E ), (60) 



where 



(0), , 16\/27r 3 a s (/it 

A s {M - - 



G F V tb V t 



ts 



Cf b ]g{Hb) H C^fg(fi b ) H -C ( 7 °Jg(n b 
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(61) 



BR(b -> s 7) for the SM 



m t rns. 


BR(b -> s 7 )[10- 4 ] 


174.3 GeV 1.6 GeV 0.3003 

171.4 GeV 1.6 GeV 0.3003 


3.32 ±0.19 
3.27 ±0.19 


171.4 GeV 0.2345 GeV 0.22 

174.3 GeV 1.6 GeV 0.22 

171.4 GeV 1.6 GeV 0.22 


3.71 ±0.20 
3.61 ±0.20 
3.59 ±0.20 



Table 3: Input values for the evaluation of BR(b — > s 7 ). For these cases, we have used a s {mz) = 
0.1185. The value of E~ = 0.2345 GeV corresponds to 5 = 0.9. The error quoted is just the 
parametric error. 

and X~ (/ife) is the analogous contribution from the primed operators (of flipped chi- 
rality L *->■ R ) of Eq. fl52l . The bremsstrahlung function B(Eq), also receives gluino 
contributions [23], which are added to the Wilson coefficients at fit, 

Cj 0)T °V) = cP 8M M + CQ*^M+xF>M,i = 7,8, (62) 

the contribution from BSM processes to the rest of the Wilson coefficients, i = 1,..,6, 
is negligible. We will denote by Kj eS g the contribution to BR(b — > s 7) containing the 
Wilson operators from gluinos: 

kl s - g = P(E )- g + N(E ). (63) 

As stated before we have implemented the calculation of BR(b — > s 7 ) as introduced in 
|22j . and adopted in [23] for the contribution of beyond SM effects. The main difference 
in the way BR(b — > sj) is calculated in the SM as introduced in [22], with respect 
to previous analyses is that m v ° le / 'm^ oZe is replaced by m^ s (m c )/ml° le in the element 
(X s j\(sc)v-A(cb)v-A\b)- In addition, the running of the charm and top contributions from 
the EW scale to the decay scale /x& is different. In Tabled we compare the values obtained 
in our calculation for the different choices of rn v ° le j 'm\ ole = 0.22 and m^ s (m c )/ml° le = 
0.3003. Lines 1 to 5 correspond to the calculation of BR(b — > s 7) implemented here 
following [22]. Line 5 is the SM value that we obtain in the decoupling limit. Line 4 is 
what we obtain when using the top mass assumed in [22] , and is in perfect agreement with 
the value obtained there (Eq. 4.14 of that reference). Line 3 corresponds to the SM when 
5 = 0.9 shown for comparison to other works. Lines 1 and 2 assume m v ° lm p h ole = 0.22 
for two values of m t . The remaining parameters for the calculation of BR(b -> s 7) that 

we use and that were not specified before are BR(B — > X c e^) exp = 0.1059, 
and /it = 5 GeV. 



0.97 
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B Interactions of fermion, sfermions and sparticles 



B.l Quarks-squarks-gluino interactions 

The interaction between quarks-squarks and gluinos is described by the Lagrangian: 



- C - - 



a/3 



-q'i'PLgaq'm + q^VrMu + h - c 



(64) 



in the electroweak-color basis, where Vr : l = (1± 7s)/2, a is the gluino color index, a,/3 
are the quark-squark color indices and % is the generation index. In the gauge basis, the 
couplings at the quark-squark-gluino vertex are given by, for incoming gluinos to sfermion- 
fermion, C-j^ 1 (y), = —igsV^T = 111 this basis the effective mass Lagrangian 

that we construct from the soft Lagrangian of Eq. (JT1) is 



£ 



cff 



2 V- 

q'JiJ 



q£ 



(65) 



where 



(M 2 f ) 



fj% 3 



ML 
M 2 LR 



M 2 RR 



+ (M))y + D 



U L,R 



C os2pMl(Tf-Q fLR sm 2 6 w ) 



~( a fij v f + /^ tanP i3( M f)ij] 
(Mj^ + (Mj) l3 + D R 

j 1, f = d 
I -1, f = u. 



P 



(66) 



In Eq. (1661) i,j = 1,2,3 are the family indices, and (iW/)^- are the non-diagonalized 
fermion mass matrices. We recall that in the eigenmass formalism all the interactions 
are computed with the particles that are mass eigenstates. Hence the effective soft mass 
matrix of Eq. (1651) . with the corresponding family elements Eq. (1661) . needs to be rotated 
where the fermions are mass-eigenstates, i.e. to the so-called super CKM (SCKM) basis. 
For the diagonalization of fermions we adopt the convention 



{u L , R )i = {KlR)ij{ u 'L,R)j, {d L ,R)i = {V^ijid'L^j, -> V; 
and thus the sfermions must be rotated similarly: 

(u L ,R)i = (V^Rhiu'L^j, (d~L,R)i = {Vt] R )ij{d' L , 

Then the soft mass matrix Eq. (16"6"|) becomes 



CKM 



R)j- 



(67) 



(68) 



(Mf KM ) 2 



M sckm2 ll M SCKM t^ 
M sckm2 lr M SCKM RR j 



{V[ T M^V[*)ij 
/t 



•((VI 



+ Ml + D{ 



((Vi T a f V^) ijVf + ^t^pM fi ) 



rf* 



f)ij 



+ M 2 fi + D f R 



(69) 
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where _D{ r remain diagonal, and Mf is the diagonal matrix of the / type fermions. The 
sources of flavour violation then come from 



SCKM 2 



{ a 'f 



(m ; 



SCKM\ 

' ij 



SCKM^ 



fa 



1.1 



{vi T Mivt: 

-(Vfa f v£ 



(70) 



which are not generically diagonal in the SCKM basis. For the case of 1 family the matrix 
{Mj u )fj is the usual 2x2 soft mass matrix whose diagonalization is obtained through 



}li f. 



R 



M SCKMZ LL M SCKM ^ LR 

M SCKM2 M SCKM2 ' 



ft 

h 



/Lj f 



R 



71 



LR RR J L 

where then the sfermion mass eigenstates, ^ are defined by 



J2 



K 



ft 

h 





= K T 


' fi ' 


. fx . 







(71) 



(72) 



Analogously for three families, we have 



11 



2\t 



(M 



<12 



'12 



(M 



'13 



(M% 3 



u 



' u L ' 




d L 


UR 




d R 




, d' = 




cr 




sr 


h 




h 


in 







(73) 



Hence the effective soft mass matrix (Ai 2 ^)ij is a 6 x 6 matrix, and the electroweak eigen- 
states and mass eigenstates are related through 

/' = (K f ) T f, (74) 

/' and / here are column vectors. The primed states are the eigenstates of the Lagrangian 
in the gauge basis and the unprimed in the mass eigenbasis. Using the above definition of 
Eq. (1721 and Eq. (1741) the couplings of the mass eigenstates involved in gaugino/ sfermion 
vertex in terms of the gauge eigenstates are given as follows 



C 9f*{fL)k 



C gf;(f R )k 



K j(2k) C gf : 



(2k) 



(75) 



We are working in the basis where Y u is diagonal at the GUT scale. At the EW scale 
two loop corrections will make it non-diagonal but due to the smallness of the expansion 
parameter e u , we are neglecting this contribution. Then Vckm — Vl, contrary to the 
convention used in many other papers (e.g. [141 [T9] ) where is taken to be diagonal 
at the EW scale. Hence the couplings of fermions and sfermions to other supersym- 
metric particles must be defined accordingly. These couplings are given in the next two 
appendices. 



38 



B.2 Quarks-squarks-neutralino interactions 



The interaction between quark-squark and neutralinos it is described by the Lagrangian, 
in the mass eigenstate basis, by: 



— djdk 



"P -I- "P 



dj > Xi 



y° \° 



(76) 



The mass Lagrangian is given by C = —\%l) NN^M.^ N*N T {%l)°) T + h.c, where the gauge 
eigenstates are ip° = [—ib, —iw, hd, h u ], for /i 12 = h±u- The mass matrix of neutralinos 
is diagonalized by M x o = JSftM.jpN* , and hence the mass eigenstates, Xii are given by 

X? = $°N)i. 

The couplings of incoming neutralinos to sfermion-fermion in the gauge basis, Eq. ( 1761) . 
can be identified from those of Eq. D20 in [39] as C,i fi (f r \ = a,j, f, (f ,\ , and 

Cxoti ici\ = br ii t f i\ , for p — L,R. Then the coupling at each vertex in the basis of 

mass eigenstates can be written as —i(C UR VL + C Ul Vr) and —iC~ x {Cd R VL + Cd L Vn\ for 
up and down quarks respectively. Cf R is any coupling involving a right-handed fermion 
state, for example C^ j^j R ^ m . Then the couplings of the mass eigenstates involved in 
neutralino or chargino/ sfermion loops in terms of the gauge eigenstates are given as 



9i 



^ijk 




r\u _ 
'-'ijk ~ 


= C X?*u){u L ) k = 


(~\d 

^ijk 


= x°M^ R )k = 


U ijk 


= C xTdMl) k = 



r d 91 



N*iKj(2k)~^yd R + ^3i^j{2k-l)^kk 



(77) 



B.3 Quarks-squarks-chargino interactions 

The quarks-squarks-chargino interactions are given by 



= ujd k 

+ *te 

+ djUk 

+ fc 



c 'x-^(4) fc ' PL + a ' +CJ ^ t ^ 



x+%(4i)fc 



C — c j / t n ^fl + C[ +c - , t , Vl 

Xt rfj«)fc Xj 



^ j ) X« 



Xi 5 Xi 



(78) 



The mass Lagrangian is given by C = — \ ifj + VV + M^+U*U T (ip ) T + if) M^ + (ip + ) T 
where ip + = [— iw + , h+], if>~ = [— iw~ , h~f], h\£ = h d , u - The mass matrix is diagonalized 
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by M^± = V + M.^+U* = U + A4^ + V* , hence the rotation to chargino mass eigenstates is 

x + = 4>l)v, x = (fc ^")^- 

The couplings of incoming charginos to sfermion-fermion in the gauge basis can also 
be identified to those of Eq. 20 of [39] . Using the transformations described in this section 
to go to the eigenmass basis we have 

3 





= C x+d*(u L ) k 


= ^E[4" 

n=l 


-l)(KlKM)nA;( — ^kk ) 




D ijk ~ 


= C x7 C d*(u R ) k 


n=l 


-i)^YcKydnk9 + U 2 , 


Kj(2n)( ^nn)(^CKM)nfe 


rtd 

D ijk " 




3 

= [ V 2i K j(2n) 
n=l 




gViiKJ( 2n _^ (V CKU )kn 


\ d , 


_ rj 


3 

n=l 


-l)( _ ^l)(^CKM)fen 





where it is useful to remember than in this notation the components Kj( 2n -i) are related 
to left-handed transformations and Kj 2n with right-handed ones. 
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